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ABSTRACT 

Context. This is the third in a series of papers that deal with angular momentum transport by internal gravity waves. We concentrate 
on the waves excited by core convection in a 3 Mq, Pop I main sequence star. 

Aims. Here, we want to examine the role of the Coriolis acceleration in the equations of motion that describe the behavior of waves 
and to evaluate its impact on angular momentum transport. 

Methods. We use the so-called traditional approximation of geophysics, which allows variable separation in radial and horizontal 
components. In the presence of rotation, the horizontal structure is described by Hough functions instead of spherical harmonics. 
Results. The Coriolis acceleration has two main effects on waves. It transforms pure gravity waves into gravito-inertial waves that 
have a larger amplitude closer to the equator, and it introduces new waves whose restoring force is mainly the conservation of vorticity. 
Conclusions. Taking the Coriolis acceleration into account changes the subtle balance between prograde and retrograde waves in non- 
rotating stars. It also introduces new types of waves that are either purely prograde or retrograde. We show in this paper where the 
local deposition of angular momentum by such waves is important. 
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1. Introduction 

Internal gravity waves (IGWs) have received growing atten- 
tion in the past 15 years as a source of angular momen- 
tum redistribution, after the suggestion by various authors that 
they could be responsible for the Sun's quasi-flat rotation pro- 
file (Schatzman 1993; Zahn, Talon & Matias 1997; Kumar & 
Quataert 1997). Charbonnel & Talon (2005) convincingly show 
that IGWs can play a major role in carrying angular momentum 
from the rapidly rotating core, left after the Sun's original con- 
traction, to the surface convection zone that is continuously spun 
down by magnetic braking. In self-consistent evolutionary mod- 
els that also takes into account the transport of angular momen- 
tum and chemical species by meridional circulation and shear 
turbulence, they also showed that the associated destruction of 
lithium is compatible with the destruction measured in the Sun. 
In the case of F dwarfs. Talon & Charbonnel (2003) show that 
the combination of surface breaking for stars with an effective 
temperature T^g < 6900 K and the appearance of IGWs at a 
temperature Teff < 6600 K when the surface convection zone is 
thick enough can together explain the existence of the so-called 
lithium dip discovered by Wallerstein et al. (1965). 

Now the question arises as to whether IGWs generated by 
a convection core also play a role in the evolution of massive 
star£]. Indeed, rotational mixing is now admittedly a major in- 

Send offprint requests to: Suzanne Talon 

' By massive stars we mean objects for which central hydrogen burn- 
ing occurs in a convection core, i.e., stars with Mi„i > 2.2 Mq at solar 
metallicity. 



gredient in such objects (see e.g. Maeder & Meynet 2000; Heger, 
Woosley & Langer 2000), and the existence of another major 
source of angular momentum redistribution within the star must 
be examined. 

However, as a general rule, massive stars are fast rotators. 
Here, fast must be understood in comparison with the wave 
frequencies that dominate angular momentum redistribution in 
stars; that is, cr 1 ^uHz (Talon & Charbonnel 2005). Our goal 
in this paper is to establish how wave properties are altered by 
rotation in preparation for full calculations of angular momen- 
tum redistribution in massive stars. §|2]is devoted to examining 
the (horizontal) eigenfunctions in the presence of rotation, while 
§[3]presents the results on angular momentum transport itself. 

2. Eigenfunctions in rotating stars 

2.1. The traditional approximation 

In the study of internal waves in the presence of rotation, two 
effects should be considered: the star's distortion caused by the 
centrifugal force, and the modification to the momentum equa- 
tion by the Coriolis acceleration. In 1-D modeling, the stellar 
distortion can be accounted for by calculating an average cen- 
trifugal force that locally reduces gravity. This effect is quite 
important at the stellar surface where the ratio r^Q^/GMr is at 
its maximum, and it modifies the location of the star in the HR 
diagram (see e.g. the ZAMS position of stars of various veloci- 
ties in Fig. 5 of Talon et al. 1997). It does however have a mi- 
nor impact on stellar evolution, since it modifies the physical 
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Fig. 1. Eigenvalues of Laplace's tidal equation when rotation is 
slowly increased. Illustrated modes correspond to the £ - 6 non- 
rotating case. 

conditions near the stellar core only slightly (see e.g. Maeder 
& Meynet 2000). In this work, we will neglect this term. The 
Coriolis acceleration, on the other hand, has a direct effect on 
the oscillations and modifies the momentum equations for the 
displacement ^, which become 



■ (T^r - 2/crQ sin 
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where cr is the wave frequency, and other symbols have their 
usual meaning. Here, we have neglected variations in the gravi- 
tational potential (<!)' = 0, Cowling 1941), which is particularly 
justified for gravity waves that have mostly horizontal displace- 
ments. This equation was given in Talon (1997), and it differs 
from the one used by other authors (see below) only by the 
derivative in Q.. 

To proceed further, we make a simplification in our system 
of equations ([1] |2l [3]) and neglect the horizontal component of 
the rotation vector = Q. sin 0. In the case of Eq. ([TJ, this is 
justified by comparing this term with the vertical gradient of the 
pressure perturbation. In Eq. (O, we invoke the dispersion rela- 
tion of gravito-inertial waves (see e.g. Unno et al. 1989) 



N^kl + (2Q ■ ky 
1? 



(4) 



which in the case cr^, (2Q)^ <K A^^ (the condition for the approx- 
imation to be valid), yields 



» kl. 



(5) 



This is similar to the result obtained for pure gravity waves. 
Since these low-frequency oscillations are quasi-incompressible, 
Eq. (|5]) is equivalent to 




Fig. 2. Top: Spherical harmonic Yf with ^ = 5, m - +2 Bottom: 

Corresponding Hough function Qkms'""^ with k = -1-3, m = +2 
and V X 2. 

whence the neglect of the second term (oc of Eq. (O. This is 
known as the traditional approximation, which is well known in 
geophysics (see e.g. Eckart 1960), and it was used for the first 
time in an astrophysical context by Berthomieu et al. (1978). 

In this case, wave eigenfunctions can be separated into ra- 
dial and horizontal components. For the radial displacement, for 
example, one get^ 

^,(r,0,0,f)=^rW0(0)e""V'^'. (7) 

The longitudinal eigenfunction is no longer given by an associate 
Legendre polynomial, as in the non rotating limit, but obeys in- 
stead a new differential eigenvalue problem of the form 



Xv [0(0)] = -A© (0), 
where 



dx \1 



(8) 



(9) 



1 



-I- mv- 



1 + v^x^ 



1 



In this equation, we use the rotation parameter 
2Q 



(10) 



(6) 



- Note that here we are using the convention of negative m for pro- 
grade modes, contrary to what has been used in previous papers. 
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and define x - cos 6. The azimuthal order takes integer values 
m - 0,+\,+2, . . . This equation, known as Laplace 's tidal equa- 
tion, was originally derived in the context of terrestrial tides. 
It forms a Sturm-Liouville problem, ensuring the existence of 
a base of orthogonal eigenfunctions. These were described in 
the geophysical context by Hough (1898), and the first numer- 
ical solution was calculated by Longuet-Higgins (1968). More 
recently, they have been examined in the context of neutron 
stars (Bildsten, Ushomirsky & Cutler 1996) and in studies of 
low-frequency oscillations in rotating stars (Lee & Saio 1997; 
Daszynska-Daszkiewicz, Dziembowski, & Pamyatnykh 2007 
and references therein). We wrote a relaxation code to solve this 
equation numerically (see Press et al. 1992). A similar procedure 
is also described in Lee & Saio (1997). 

The Hough function @km (v, cos 9) is associated with the 
eigenvalue Kkmiy)- Following Lee & Saio (1997), the order 
k = 0, + 1, ±2, . . ., which does not explicitly appear in Eq. (|9]l, is 
such that Ak+i,m (v) > Afa„ (v). In the v = limit, Eq. (8) is iden- 
tical to Legendre's equation with k - {-\m\ and Afa„ - £({ + I). 
Slowly increasing v, we first lift the degeneracy in A (see Fig.!!). 
Increasing rotation confines eigenmodes toward the stellar equa- 
tor (see Fig.|2]i. 

In the case of more rapid rotation (i.e., cr < 2Q), mode or- 
ganization becomes much less trivial. First, new solutions ap- 
pear when V > 1 and have negative eigenvalues; thus they are 
only evanescent. However, as v increases, so do their eigenval- 
ues, which become positive for the retrograde modes (and so the 
waves become propagative) when v - (m - k) (m - k - I) /m 
(cf. Lee & Saio 1997). Second, a large fraction of the modes 
that exist in the absence of rotation see their eigenvalue grow 
tremendously with rotation (see Fig. [3]l- Townsend (2003) has 
introduced a new classification for modes based on the behavior 
of their eigenvalues in rapid rotation: 



- Gmvito-inertial waves: they are similar to normal gravity 
waves, but modified by the Coriolis acceleration. They are 
assigned indexes s - 1,2,3,..., such that As+i_,„(v) > 
A^m (v). 

- Rossby waves: these are purely retrograde waves, which ex- 
ist only in the case of rapid rotation. They arise from the 
conservation of specific vorticity, combined with the effect 
of curvature. They also have indexes s - 1, 2, 3, . . . but or- 
dered such that As+i_,„(v) < A,,,,, (v). The s - I modes are 
retrograde sectoral waves. 

- Yanai wave^ they behave like a mixture of gravity and 
Rossby waves. The m < modes exist in the absence of ro- 
tation. The m > modes appear when v - m + 1 with small 
eigenvalues, and their horizontal eigenfunctions are then ex- 
actly 0(v = OT H- 1) - P'"j^i- When they appear and have 
small eigenvalues, they behave mostly like Rossby waves; 
m < waves and m > waves with large eigenvalues behave 
rather like gravity waves. We assign to them index 5 = 0. 

- Kelvin waves: like Rossby waves, they arise from the conser- 
vation of specific vorticity combined this time with the strat- 
ification of the medium. They are purely prograde waves, 
whose characteristics change little with rotation. This occurs 
because their displacement is very small in the 6 direction. 
They correspond to the m < Q sectoral modes. We assign to 
them index 5 = - 1 . 



C 




Fig. 3. Eigenvalues A of Laplace's tidal equation for rapid 
rotation (top) and equivalent eigenvalues A for the horizon- 
tal Laplacian of Hough functions (bottom). Mode classifica- 
tion: Black: Gravito-inertial waves Blue: Yanai waves Brown: 
Kelvin waves Green: Rossby waves. Continuous lines corre- 
spond to numerical solutions and dashed lines to asymptotic so- 
lutions (see § 12. 3t . 

Without rotation, the correspondence between various indexes 
follows 

s - £ - m + \ for m > (11) 
and 

s ^( + m-l form < (12) 
and is illustrated in Fig.|4] 

2.2. Horizontal displacement 

Pressure (P') and density (p') perturbations have a horizontal 
structure identical to that of the radial displacement (Eq.|7J. This 
is, however, not the case for the two components of the horizon- 
tal displacement, which take the form 

1 P'(r) 
rcr^ p 

i P'(r) 

^^(r,0,0,f) = - 
rcr^ p 

with 

•7/1 (y, cos 0) 



'H^v,cos0)e""*e"^' 
'HHv,cos9)e""'^e''" 

1 



(13) 
(14) 



^ They have been named after their discoverer in the Earth's atmo- 
sphere, Yanai & Maruyama (1966). 



(1 - v^x^) ^Jl-x^ 

— 1 — X h mvx 

^ ' dx 



@sm (V, COS 9) 



(15) 
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Fig. 4. Correspondence of index s with index i in the absence of rotation. Modes in the dashed box (retrograde Yanai waves) do not 
exist in the absence of rotation. 




-1 
1 

C.ri 
-0.5 





A : 






^ — \ — \ — \ — \ — — \ — \ — \ — \ — 


— 1 — 1 — 1 — 1 — — 1 — 1 — 1 — 1 — ^ 













COSI 



Fig. 5. Hough functions 0^,,,, "Kf,,,, and ^fm- Normalization is given by d0/d cos = 1 at cos 0-0. Left: s - 4-, m - 2, and v = 
({ = 5). In this case, 0,„,„ "Kf^, and "Hf,^ ai-e given by P"/, dP"//de, and mP™ / sin 6 respectively. Right: 5 = 4, m = 2, and v ^ 2. 



•Ht(v,cos0) = 



(1 - v2x2) vr^ 

—vx(l — x^) h m 



V, COS 6) 



(16) 



where again x - cos 0. An example of these functions is shown 
in Fig.|5] Note that ©j,,, has two more zeros compared to the non- 
rotating case. Furthermore, 'H^,„ and '7/f„, are no longer given 
by simple expressions and have zeros that do not coincide with 
those of 05„. 



2.3. Asymptotic solution 

When considering the impact of IGWs on the distribution of 
angular momentum within a star, low-frequency waves (with 
cr XI I jjHz) play a dominant role (this is discussed at length 



in Talon & Charbonnel 2005). In the case of massive stars that 
are generally rapid rotators (with Q » 20 fjHz), we are within the 
limit V » 1, where asymptotic solutions to Laplace's equation 
exist (Townsend 2003). 

To get these solutions, Townsend defines 



(cos 61) = sin fl'K* (cos 61) 
(cos 0) = - sind'H'^ (cos 6) . 



(17) 
(18) 



Then, eliminating 0, Laplace's equation is equivalent to two 
coupled first-order equations 



[l-x^) — -mvx = (A^-l)0 



[l-x^) — + mvx = [a (l - x^) - 0. 



(19) 
(20) 
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Further simplifications to these equations can be made. For 
gravito-inertial waves and Yanai waves, A » m^, while Rossby 
waves have K (see Fig. O. When v » 1, these waves 

are very confined to the equator, which allows the assumption 
^1 - x^) =s 1. For Kelvin waves, one has A =s m^, leading to 

|0| » |0|. Townsend (2003) then finds approximate solutions 
to the Hough functions in the form of Hermite polynomials. We 
refer the reader to his paper for their derivation. Let us note here 
that, except for Rossby waves that are less confined to the equa- 
tor, a direct comparison with the numerical solution for the com- 
plete version of Laplace's equation shows that the approxima- 
tion is valid for v > 2. For Rossby waves, this approximation 
becomes valid at much higher values on the order of v > 20 (see 
Fig.E). 

In this asymptotic limit, eigenvalues take the form 



A„„ (v) 



-v(2s + l) + 
2 ' 1 



(ls+ Yf -Mmv-rn^) 



for gravito-inertial waves (s - 1,2,...), 

t^sm (v) = (v - mf- 

for Yanai waves (s - 0), 



l^sm (y) = m 

2mv + 1 

for Kelvin waves (s - 



-1), and 



A^m (v) = 



-v(2s H- 1)- - -y/y2 (2s -H 1)^ - 4 (mv - m^) 



(21) 



(22) 



(23) 



(24) 



for Rossby waves (s - 1,2,...). These asymptotic values are 
shown in Fig.|3]and are given for v > 2 for gravito-inertial waves, 
Yanai waves, and Kelvin waves and for v > 20 for Rossby waves. 
As V increases, so does the agreement between the exact and 
asymptotic value of A. 

For gravito-inertial, Yanai and Rossby waves, eigenfunctions 
are given by 

(VAv)'^^ 



K-nP- 

(m \ I I m 



-iU/.v+i(0 



(0 = m 



(VAv) 



1/2 



A — irP- 



(25) 



(26) 



(27) 



fVA 

s\ + 1 

1 m 



^.-1 (0 + 2 



with ( = cos 6 ( VAv) ^ and where H, is the Hermite polyno- 
mial of order s. For the Kelvin waves, we get 



rvA 



- 1 



-r/2 



0(7?) = e 



0(/7) = - 



-r/2 



1 



(-mvy^ 2mv + 1 

2 



1-/2 



Of) = —ffi 



+ l\ rie 



rje 



\ 2mv + 1 

with 77 = cos 0(-OTv)''^ (Townsend 2003). 



(28) 
(29) 

(30) 



2.4. Radial dependency 

The formalism used to treat the wave radial dependency is not 
modified compared to the case of the non-rotating star. In the 
anelastic approximation cr^ <s: N^, the adiabatic radial displace- 
ment ^r{r) obeys 



dV IN^ ,\A 

-t:t+ -T-iK-A^o 



(31) 



with iA('-) = ^/pr%■ (r), A^^ = A^^ + ^2 ^ g^y^^ _ + 

and where we neglected the righthand side term (cf. Press 1981). 

The radial wavenumber is still 



- 1 



A 
3- 



(32) 



However, in contrast to the non-rotating case, the horizontal 
wave numbeiQ is no longer given by - A/r^. When horizontal 
gradients are required, an equivalent value has to be taken that 
we define here as 



r2v20,,„(y)p)^ 
(|0.m(v)l')^ 



(33) 



(see Fig.|3]bottom), where (. . .}g denotes latitudinal averages. In 
the absence of rotation, we recover Asm = A „„. We may define an 
'equivalent' horizontal wave number that is now given by k^^ = 
A/r^, and we have kj » kj^ as in the non-rotating case (except 
when A — > 0). 

We now apply the WKB method to solve Eq. ( ISTT i. For 
a wave propagating towards the surface, we write iff{r) = 
C(r)e-'' /'''■'*' and get 



^Ar)oc\k,r"'p-^'^r-^e-<I'^'''': 



(34) 



Assuming that wave damping is dominated by heat diffusion and 
that this term remains small, the radial displacement becomes 

^,(r)e-^^'^'^ with 



T(r) 



A3/2 KNNj I 



(jA Y^2 



(35) 



(Zahnetal. 1997). 



3. Angular momentum transport 

3.1. Excitation 

The treatment of wave excitation is more complex in the case 
of the rotating star. Indeed, the traditional approximation, which 
requires cr^,(2Q)^ <K A^^, is no longer valid in the convection 
zone; thus, formalisms that convolve the wave eigenfunctions 
with entropy perturbations and Reynolds stresses in the convec- 
tion zone (such as those by Goldreich, Murray & Kumar 1994, 
or Balmforth 1992) require major transformations (Belkacem et 
al. in prep.) For that reason, we resort to adapting the Garcia 
Lopez & Spruit (1991) formalism to our problem. 

In that model, waves are excited at the boundary of the ra- 
diative zone by the interface deformation caused by convective 



* In the absence of rotation, Legendre polynomials follow 
r^V^© = -A0, 

but this is no longer the case when rotation is added (see Eq.[9ll. 
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Fig. 6. Energy flux in IGWs at the base of radiative envelope in a 
3 Mq main sequence star with Xc - 0.5. Black: Gravito-inertial 
waves Blue: Yanai waves Brown: Kelvin waves Green: Rossby 
waves. Continuous lines correspond to numerical solutions and 
dashed lines, to asymptotic solutions (see § 12. 3b . 



eddies. The eddies mean energy density is transfered to waves of 
the same frequency and with ki, < ^eddy according to 



(36) 



(Garcia Lopez & Spruit 1991), where (...)/, denotes horizon- 
tal averages. The factor kh/ks^dy reflects the reduced efliciency 
in the production of small ki, waves caused by the stochastic na- 
ture of turbulent eddies. Assuming a turbulent Kolmogorov spec- 
trum, we have Leddy °^ '^Iddy < Lc with Lc the mixing 
length. The spectral distribution becomes 



k^ 

"eddy 



Vddy 



(cT)^ki\ — 

CTr 



(37) 



(38) 



with cr > (Tc, kc = 2n/Lc and where v^ and cr^. are the char- 
acteristic convective velocities and frequencies, respectively. We 
get 



(39) 



To proceed further, we assume that the horizontal dependence of 
IGWs at the radiative boundary is given by Hough function^. 
Using the 'equivalent' eigenvalue given by Eq. ( [33] l, one then 
has 



(40) 



^ This is formally wrong since at this boundary we have — > 0, 
which contradicts the traditional approximation. However, as A'- rises 
rapidly away from the boundary, and considering the inescapable pres- 
ence of a slight amount of overshooting, the approximation is accept- 
able. 




Fig. 7. Ratio m' - for the gravito-inertial waves when the 

Coriolis acceleration is neglected in the calculation of the angu- 
lar momentum flux. Continuous lines correspond to numerical 
solutions and dashed lines to asymptotic solutions (see § |2.3l l. 



An alternative definition of the equivalent eigenvalue A - 

similar results. Following 
use a projection of the 



leads to very 
Mathis (2005), one could also 
Hough functions onto spherical harmonics for which horizontal 
wavenumbers are well-defined. However, in the case of fast ro- 
tation numerical problems are encountered with this last method 
because Legendre polynomials are not a good basis for project- 
ing Hermite polynomials. 

3.2. Energy flux 

Angular momentum transport by IGWs is dominated by the low- 
frequency waves. Here, we will thus adopt cr^ <s: and use 
k^. = (N/a-)^ A/r^. The kinetic energy flux per unit frequency is 
given by 



where flie radial group velocity Vg is given by 

dcr cr 
" d^ ~ "k/ 



(41) 



(42) 



In the case of main sequence massive stars, we have k,- < 0, 
corresponding to waves traveling from the convective core to the 
surface. Using Eq. ( |40] | to get the energy of a given mode, we 
obtain 



Te (o", s, m) 



2 2 



(43) 



where A^^, which would formally be at the interface, is taken a 
fraction of a pressure scale height into the radiative region. 

These fluxes have been computed for a 3 Mq main-sequence, 
population I star calculated with the Geneva stellar evolution 
code and for a rotation velocity of Q = 20yuHz (Fig. This 
value, typical of massive stars, corresponds to ~ 40% of the sur- 
face critical velocity. 

In the case of Rossby and retrograde Yanai waves, one can 
see that they appear with a significant energy flux. This is caused 
by the divergence of the group velocity when A„„ is close to 
(see Eq. l42b . This situation is unphysical; however, it will be 
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cosff 

Fig. 8. Left: Hough functions 0s„, and ^t,„ for 5 = 4, m = 2 
angular momentum flux. Normalization is given by d0/d cos 9 - 
the rotating case, the product sind&'N'^ changes signs locally. 
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= 5), V = 0, and their product sinOQ'H^, which is used for the 
1 at cos 6-0. Right: Equivalent Hough functions for v ^ 2. In 




5 10 15 20 20 40 60 80 



Fig. 9. Ratio m' = -f ^ for various modes. Top left: Gravito-inertial waves corresponding to ^ < 6 Top right: Yanai waves of 
order m - -5, . . . , +3. Bottom left: Kelvin waves of order m > -6. Bottom right: Rossby waves of order s - 1,2, and m < 4. 
Continuous lines correspond to numerical solutions and dashed lines to asymptotic solutions (see § l2.3l l. 



compensated for by a term in A in the calculation of the angular 
momentum flux (Eq.l48]l. Thus we do not need to correct this in 
the present treatment. 



3.3. Angular momentum flux 

We must now convert the energy flux into an angular momentum 
flux. Following Zahn et al. (1997), we write 



(44) 
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with the real part of the velocity field v = icr^ 
V, = A(r) sin lo-f - I Ar'k,. + mM (v, cos 6) e'^^'^'^ 



-Jdr'kr 



rkr 

— -A(r) cos I erf 
A 



rk. 



- dr'kr + m^j fi" (v,cosff)e 
■A(r) sin |crf - dr'kr + 'H'* (v, cos 0) e"^ 



given by the WKB method for Vr, and with help of the continuity 
equation in the Boussinesq approximation for v^ and v^. The 
energy flux for a given mode becomes 

^.=-^.(K>.^f [(i«f>/(i-«f>J)r-'" 

This equation is equivalent to the one obtained by Mathis (2005). 

In the absence of rotation, the angular momentum flux is ex- 
pressed as (Zahn et al. 1997) 

Tj = (prsin6»v,.v^)^^ 

which becomes 

l2 



(45) 



(46) 



for a given mode. When this simple formulation is applied to 
gravito-inertial waves, one finds that the sign of angular momen- 
tum transport varies with v (see Fig.|7]l. As the rotation param- 
eter V increases, the value m' = -f diminishes linearly and 
changes sign for v = 2m. This would imply that a retrograde 
wave (m > 0) would, for fast enough rotation, carry a posi- 
tive angular momentum flux (m' < 0). This can be understood 
with the aid of Fig. |8] In the absence of rotation, - F" and 
"W* = mP"y sin 6*; their product always has the sign of m. This 
is no longer the case when the Coriolis acceleration is included. 
The integral ^sin00'K'*^^ may thus change sign. 

The solution to this paradox lies in the fact that, in a rotating 
system, the angular momentum flux is actually given by 



fj = (jjr sin Vr (v^ + 2Q cos 



(47) 



(Jones 1967; Bretherton 1969). The first term in this equa- 
tion corresponds to the angular momentum flux across an 
Eulerian surface, and the second term to a flux associ- 
ated with a Lagrangian contribution to angular momentum 
(Bretherton 1969). In that case, the actual angular momentum 
flux thus becomes 



AHr) 



rk 

pr ((sin e "K'*)^ - V (sin cos 6» -H") J — e"^'''- 



A 



for a given mode. The ratio between the kinetic energy and the 
angular momentum flux is thus 



Ty 2A (sin ^ 0'^'^)g-^ (sin f COS 61 0'H'')^ 



rE 



where we use the approximation r^kf. » A^. In the absence of 
rotation, we get 



2m 
a 



(49) 
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Fig. 10. Local angular momentum luminosity normalized to 1 at 
the base of the radiative zone (top) and local deposition of an- 
gular momentum (bottom) for waves with cr = 3 /iHz in a 3 
main sequence star (Xc - 0.5). Black: Gravito-inertial waves 
Blue: Yanai waves Brown: Kelvin waves Green: Rossby waves. 
Continuous lines correspond to numerical solutions and dashed 
lines to asymptotic solutions (see § |2.3l l. 



We define an equivalent azimuthal number m' such that 

o-^r^ (sin00'K'*) -v(sin6'cos00'K'') 

^ = -ttt;^ = ^ ; — (50) 

2Te /i^"|2\ • ■ 



(l^*l')« 



which has been evaluated numerically and is shown in Fig.|9l 

In the case of gravito-inertial waves, as v increases, the value 
m! rapidly converges towards m/3, for all orders s. This implies 
that the symmetry that exists between prograde and retrogrades 
waves is conserved in the rotating cas^ Yanai waves behave 
similarly, but with a slower convergence rate. 

In the case of the Kelvin waves, m' varies only slightly with 
rotation and remains close to m. Their angular momentum flux 
is always positive. Rossby waves appear with m' ^ 0, and its 
value rises slowly and tends towards m with the increase in the 
rotation rate. Their angular momentum flux is always negative. 



(48) 3.4. Angular momentum deposition 



The angular momentum distribution within the star evolves un- 
der the effect of the damping of IGWs. When traveling inwards, 
each wave deposits its angular momentum at the location where 
it is damped. In stars, the major source of damping is thermal 



Note however that excitation could be asymmetric. 
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diffusivity {cf. § |2.4| i. which is different for each wave. We de- 
fine the angular momentum luminosity £.j - Anr^Tj, which in 
the WKB approximation can be expressed as 

j:j{r)^An?Tj{rc)e-'^'-\ (51) 

Then, the local angular momentum evolves according to 



d /Stt ^ „ , 

d^ -^"'■ ) = -d; 



\ d i r 



(52) 



(Zahn et al. 1997). 

In this exploratory step, we examine where damping occurs 
for each type of wave. Let us first look at Fig.[TOl which gives the 
local angular momentum luminosity and the angular momentum 
deposition -dXy/dr of waves with a frequency cr = 3 yuHz, for 
the 3 Mq model of Fig. |6] (the rotation parameter is v = 13.3). 
Since all waves shown here have the same frequency, the loca- 
tion of the deposition is determined by the eigenvalue A (see 
Eq. [35] l. The rapid fall in the local amplitude at the convec- 
tive boundary (over 6r ^ 1%) is caused by the mean molecular 
weight gradient that has been left behind the regressing core, and 
damping increases in the outer region with the thermal diff'usiv- 
ity K. 

4. Discussion 

In this paper, we have examined the transport of angular momen- 
tum by low-frequency waves excited by core convection. The 
traditional approximation was used to evaluate the impact of the 
Coriolis acceleration on waves. We find two main effects: 

- the horizontal structure, and hence the amount of angular mo- 
mentum carried by a wave, is modified; 

- new types of waves appear 

In the case of gravito-inertial waves, the main effect of the 
Coriolis acceleration is to confine the wave towards the equa- 
tor. This reduces the total angular momentum carried by a wave 
of a given amplitude and frequency (Fig. |9]l. As rotation in- 
creases, a saturation exists in this reduction, which tends to- 
wards ;«' - -f ^ = f ■ Rotation also increases their radial 
wave number (see Eq. [32] and Fig. |3]l, and hence their damp- 
ing. These waves are thus deposited very close to the convection 
core (see Fig. [TO] i. If the asymmetry in the excitation of pro- 
grade and retrograde waves is not too strong, we expect that the 
damping of these waves could produce a shear layer oscillation 
(SLO) similar to the one present in slowly rotating stars (Talon 
& Charbonnel 2005 and references therein). 

The Yanai waves, whose main restoring force is also gravity 
for high values of the rotation parameter v = 2Q/cr, have the 
same limit value for m' . Their eigenvalues A are, however, lower 
so they are damped farther from the convection core (Fig. [TOl i 
and over a larger portion of the star. These most certainly do 
not produce a second SLO but could generate a local shear in 
the interior (corresponding to the stationary solution of Kim & 
MacGregor2001). 

The other two types of waves, namely Kelvin and Rossby 
waves, show a somewhat different behavior. For these waves, the 
main restoring force is the conservation of vorticity, combined 
with stratification in the first case and in the second, combined 
with curvature. These two types of waves have the same limit- 
ing value m' - m, and their eigenvalues A are lower than those 
of both gravito-inertial and Yanai waves. As a result, they are 
damped much closer to the stellar surface (Fig.fTOb. These waves 
could create a strong shear in that region and could thus induce 



a large amount of mixing close to the stellar surface. Complete 
dynamical simulations remain to be completed to verify these 
conjectures. 
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